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一 种 基于 黏 弹性 理论 预测 DP600 钢板 滞后 
回 弹 的 新 模型 


孙 帅 ， 骂 大 辛 
(北京 理工 大 学 材料 学 院 , 北京 100081) 


摘 要 通过 加 载 不 同 预 应 变 和 保 压 时 间 获 得 了 DP600 钢 板 的 室温 滞后 回 弹 响应 . 基于 笑 弹 性 理论 , 采用 309 MPa 恒 应 力 
室温 蠕 变 实验 , 获得 了 预测 滞后 回 弹 过 程 的 蠕 变 和 柔 量 函 数 . 基于 DP600 钢板 室温 弹性 和 塑性 加 载 响应 , 修正 了 线 黏 弹性 本 
构 方程 的 积分 下 限 , 将 柔 量 函 数 与 修正 的 应 力 率 函 数 进行 卷 积 , 获得 了 影响 滞后 回 弹 过 程 的 加 载 历史 效应 . 再 将 加 载 历史 
应 变 波 形 与 卸载 脉冲 登 加 , 获得 不 同 预 变形 量 的 滞后 回 弹 应 变 预 测 值 . 结果 表明 , 滞后 回 弹 初期 应 变 率 较 高 , 之 后 应 变 率 逐 
渐 降 低 ; 随 着 预 应 变 的 增加 , 滞后 回 弹 应 变 的 绝对 值 逐渐 递增 , 并 且 滞 后 回 弹 占 总 回 弹 比例 逐渐 增加 ; 滞后 回 弹 预测 方向 与 
当时 回 弹 方向 相同 , 与 加 载 方 向 相反 ; 相同 时 间 段 内 , 蜂 变 回 弹 应 变 比 灌 后 回 弹 应 变 明 显 降低 , 与 滞后 回 弹 实 验 相 比 , 蠕 变 
恒 应 力 与 蠕 变 回 弹 应 变 的 比值 随 预 应 变 的 变化 规律 呈现 出 良好 的 线性 关系 ; 基于 和 弹性 理论 的 修正 模型 所 得 到 的 预测 值 
接近 实验 测量 值 , 并 能 够 较 好 地 反映 实验 变化 规 得 
关键 词 DP600 钢 板 , 滞后 回 弹 , 蠕 变 回 弹 , Kelvin 模型 
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ABSTRACT DP600 steel sheet with high strength has drawn much attention in the automotive industry, but the 
shape change following forming and unloading has not been known widely. The time- dependent springback of 
DP600 steel sheet was investigated under different pre-strains by uniaxial tension. According to the viscous behav- 
iors under the elastic and plastic loading tests, the lower limit of integration in the constitutive equation of linear 
Viscoelasticity was modified and the creep compliance was gained from the creep curve at a constant stress level of 
309 MPa at room temperature. The predicted curve was acquired by using the superposition of the unloading im- 
pulse and the historical loading curve. The results reveal that the strain rates with high initial values gradually de- 
creased following unloading at room temperature. As the pre-strain went up, both the absolute anelastic strain of 
time-dependent springback and the anelastic proportion of the total springback increased. Meanwhile, the direction 
of the time-dependent springback was same as that of the initial springback and opposite to the loading direction. 
In the same springback period, the ratio of the unloading stress to the creep springback strain tended to vary more 


linearly with the pre-strain than those obtained from immediate unloading. The simulated results using the revised 
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model are in good agreement with the experimental data. 
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随 着 汽车 、 航 空 航天 种 


j 造 技术 不 断 向 高 精度 、 


轻 量化 方向 发 展 , 各 种 新 型 高 强度 钢板 和 镁 合金 


铝 合金 等 轻 合 金 板 材 的 应 用 越 来 越 广泛 , 因而 对 它 


们 的 成 形 性 特别 是 成 形 后 形状 稳定 性 的 研究 合 发 


重要 . 回 弹 是 影响 金属 塑性 
主要 因素 , 实际 上 , 很 多 金属 板材 或 管材 在 塑性 成 
形 后 会 经 历 2 种 回 弹 : 一 种 是 色 载 瞬间 时 刻 发 生 的 
瞬时 回 弹 ""; 另 一 种 是 随 着 时 间 的 
在 印 载 后 还 将 继续 发 生 微小 的 变化 , 称 其 为 “ 消 


回 弹 ”或 “时 间 依 赖 性 


人 


回 弹 的 研究 取得 了 一 


为 潜在 因素 将 直接 影响 
足够 的 重视 . 


回 弹 ” 近年 来 , 关于 金属 


emer de 而 回 弹 的 滞后 性 作 


的 滞后 回 弹 应 变 的 预测 ; 


的 预测 , 例如 弛 豫 时 间 . 


二 是 对 滞后 回 弹 特 征 时 间 


向 拉 伸 实验 有 清晰 的 应 力 - 应 变 对 应 关系 , 便 


于 研究 滞后 回 弹 规律 , 同时 也 是 研究 多 向 应 力 状态 


成 形 后 形状 稳定 性 的 最 


准 移 , 其 形状 尺 


下 应 变 啊 


应 的 重要 基础 . 然而 目前 , 对 室温 拉 伸 消 


后 回 弹 实验 现象 的 研究 
的 加 载 历 史 对 DP600 钢 


还 不 够 充分 . 为 了 揭示 不 同 
板 滞后 回 弹 的 影响 规律 , 针 


对 DP600 钢 板 进行 了 单 


些 进展 , 但 关于 滞后 回 弹 


零 部 件 的 装配 精度 , 需 引 起 


对 时 间 依 赖 性 行为 的 研究 始 于 20 世纪 40 年 代 . 
Zener% 首 次 定义 并 明确 了 金属 黏 弹 性 (viscoelastici- 
ty) 和 潍 弹 性 (anelasticity) 的 区 别 . 1972 年 , Nowick"” 


提出 , 应 将 金属 的 颖 弹性 


的 微观 机 理 
相互 作用 


与 内 耗 相 区 分 . 滞后 回 弹 
与 点 缺陷 、 线 缺陷 以 及 面 缺 陷 的 运动 和 
关 "" ,目前 , 室温 下 自然 状态 的 金属 滞 
后 回 弹 行 为 的 常见 激发 方式 为 塑性 变形 及 之 后 的 
瞬时 印 载 联合 作用 , 因而 塑性 变形 的 微观 机 理 与 灌 
后 回 弹 的 微观 机 理 有 着 密切 的 联系 . 同时 , 滞后 回 
弹 的 激发 方式 与 蠕 变 回 红 
看 成 滞后 回 弹 的 特殊 情形 . 

Wang 等 中 和 Lim 等 
昌 合 金 板材 的 单 向 拉 伸 滞后 回 弹 以 及 拉 弯 滞后 回 


弹 , 通过 调节 弯曲 半径 及 拉 应 力 等 工 


基本 相同 , 蠕 变 


回 弹 也 可 


研究 了 几 种 高 强 钢 以 及 


弯曲 角 的 滞后 回 弹 量 , 认为 在 外 载 后 2 h 内 , 主要 为 


黏 弹性 机 理 控制 . 然而 目前 


弹 的 预测 模型 还 极 少 , 看 


月 ， 专门 针对 室温 灌 后 回 


单 性 理论 的 应 


] 还 主要 针 


对 于 实验 结果 的 被 动 拟 合 上 . 


响应 , 提 


不同 的 加 载 预 应 变 , 测试 外 载 后 的 室温 滞后 回 弹 
出 了 一 种 基于 番 弹 性 理论 的 弹性 后 效 修正 


向 拉 伸 和 室温 里 变 实 验 , 采 


公式 , 重点 预测 DP600 钢 板 室温 滞后 回 弹 应 变 . 


1 实验 方法 
实验 所 用 DP600 钢 
数 ，% ) 为 : C 0.009, Si 


0.05, Fe 余 量 . 采用 WDS-10T 万 能 : 


室温 滞后 回 弹 和 蠕 变 实 


板 的 主要 化 学 成 分 (质量 分 
0.36, Mn 1.84, S 0.005, Al 
包子 试验 机 进行 
验 . 试 样 取 样 方向 平行 于 轧 


制 方向 , 平行 段 长 度 70 mm, 标 距 长 度 50 mm, 宽度 


15 mm, 厚度 1.6 mm. 引伸 计 分 辩 率 为 0.001 mm. 
拉 伸 预 应 变 为 0.06%, 0.2%, 5%, 8%, 10% 和 
15% 的 试 样 编号 分 别 为 EI-0.06, PI-0.2, PI-5, PL8， 


PI-10 和 PI-15, 其 中 首 字母 E 代 表 弹 性 加 载 ,P 代表 


塑性 加 载 ; 第 2 个 字母 大 为 I 代 表 直 接 外 载 的 滞后 


回 弹 实验 , 数字 代表 拉 伸 预 应 变 ; 若 第 2 个 字母 为 
C 代 表 里 变 实验 , 数字 代表 蜂 变 过 程 的 预 应 变 . 晰 
变 恒 应 力 取 应 力 - 应 变 曲 线 上 预 应 变 对 应 的 应 力 . 对 


于 DP600 钢板 宏观 届 


程 的 试 样 横 截面 尺寸 变 
行 分 析 . 


服 应 力 以 下 PC-0.2, PC-0.3， 
PC-0.33 和 PC-0.4 试 样 的 里 变 实验 采用 恒 载 荷 率 
500 N/s 快速 加 载 , 滞后 回 弹 实验 采用 恒 位 移 率 
1 mm/min 方式 加 载 , 各 实验 忽 


略 弹性 加 载 和 蠕 变 过 
化 , 以 恒 应 力 率 和 恒 应 力 进 


加 载 过 程 的 数据 采 


Wang 等 中 采用 室温 


蠕 变 模型 对 室温 滞后 回 弹 


进行 模拟 , 结果 表明 其 济 


论 以 蠕 变 和 Boltzmann 县 


后 回 弹 量 的 准确 性 强烈 地 
依赖 于 变形 印 载 之 后 的 残余 应 性 场 而 线 黏 弹性 理 
登 加 原理 为 基础 , 通过 卷 积 
计 入 历史 效应 对 后 期 应 变 或 应 力 的 影响 


泽 


首 且 进 一 


步 计 入 其 他 非 线性 效应 的 修正 也 更 加 灵活 "™. 因而 


吸 待 开发 基于 番 弹 性 理论 


, 能 够 反映 不 同 加 载 历史 


效应 的 滞后 回 弹 预 测 模型 . 预测 金属 滞后 回 弹 的 模 


型 主要 涉及 2 个 方面 : 


是 对 滞后 回 弹 趋 了 


平衡 时 


集 方式 采用 实验 设备 既定 


采集 频率 1.85 s 一 加 载 至 预 应 变 后 , 立即 卸载 . 当 试 


样 与 下 夹 头 完 
伸 计 所 示 位 移 作为 滞后 
蠕 变 回 弹 实验 采用 非 等 


全 脱离 时 ( 缀 载 后 20 s 之 内 ), 记录 引 


回 弹 初 始 位 移 . 滞后 回 弹 和 
时 间 间 隔 采 集 数 据 , 以 便 完 


整 记录 数据 变化 细节 . 


具体 方法 如 下 : 在 每 次 示 数 


发 生变 化 时 , 若 示 数 变化 后 稳定 , 记录 稳定 值 和 变 


化 瞬时 时 刻 ; 若 示 数 波 去 


, 取 波 动 值 的 平均 值 , 并 记 


录 其 瞬时 时 刻 . 而 滞后 


回 弹 和 里 变 实验 的 加 载 阶 


段 , 仍 按 实 验 设备 设 定 的 等 时 间 间 隔 采集 数据 . 
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2 实验 结果 与 讨论 
2.1 单 向 拉 伸 实验 

图 1 为 DP600 钢 板 单 向 拉 伸 实验 的 应 力 -应 变 
曲线 . 通过 应 力 - 应 变 曲 线 获得 的 DP600 钢 板 的 基本 
力学 性 能 如 下 : 弹性 模 量 为 182 GPa, 届 服 强度 为 
345 MPa, 抗 拉 强度 为 628 MPa, 届 强 比 为 0.549, 最 
大 应 力 塑 性 延伸 率 为 22.17%, 强度 系数 为 1062.9， 
加 工人 硬化 指数 为 0.207. 由 于 实际 生产 中 , 板 料 拉 深 
变形 基本 处 于 均匀 塑性 变形 区 , 因而 将 均匀 塑性 变 


2.2 瞬时 回 弹 

表 1 为 EI-0.06, PI-0.2, PI-5, PI-8, PI-10 和 PI-15 
试 样 直 接 秃 载 后 的 瞬时 回 弹 与 滞后 回 弹 数 据 . 可 
知 , 随 着 预 应 变 的 提高 , 卸载 应 力 Gi 增 大 , 同时 瞬时 
可 弹 应 变 Ag 也 增 大 . 定义 G 与 As 的 比值 为 沾 后 回 
弹 瞬时 模 量 E, 其 物理 含义 为 单位 瞬时 回 弹 应 变 所 
对 应 的 缀 载 应 力 ; Gi 与 滞后 回 弹 应 变 As 的 比值 为 滞 
后 回 弹 模 量 E, 其 物理 含义 为 单位 滞后 回 弹 应 变 所 
True stress-strain curve 


800 上 PI-15 
PI-10 


PI-8 
600 上 


Engineering stress-strain Curve 


400 上 上 


Stress /MPa 


200 上 


0.00 0.05 0.10 0.15 0.20 0.25 0.30 
Strain 


图 1 DP600 钢 板 单 向 拉 伸 应 力 -应 变 曲线 
Fig.1 Stress-strain curves in uniaxial tensile test for 


DP600 steel sheet 


表 1 试 样 经 直接 卸载 后 的 瞬时 
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对 应 的 务 载 应 力 : 

天 = 人 As、 (1) 

E,=0o./Ae. (2) 


将 表 1 的 EE 数据 绘 于 图 2. 可 知 , E, 随 预 应 变 & 
的 提高 而 单调 递减 , 呈现 滞后 回 弹 瞬 时 模 量 软化 现 
象 W9 . 由 于 DP600 钢 板 的 及 呈现 指数 衰减 趋势 , 因 
而 对 图 2 中 数据 , 分 别 采用 指数 与 究 函 数 进行 拟 合 ， 


E =58.41 exp(-e,/0.02832)+ 127.93 (3) 
Eb,=106.2e (4) 

2.3 滞后 回 弹 
从 表 1 中 可 见 , 对 于 弹性 加 载 实验 , 试 样 EI-0.06 
炙 载 后 , DP600 钢板 的 瞬时 回 弹 应 变 等 于 其 加 载 应 
变 , 并 且 其 滞后 回 弹 量 微小 , 低 于 引伸 计 的 分 辩 率 ， 
说 明 室 温 下 DP600 钢 板 在 弹性 应 力 秋 载 后 , 滞后 
弹 十 分 微弱 . 当 预 应 变 & 增 至 0.2% 时 , PI-0.2 试 样 的 
As 的 绝对 值 仅 为 0.006%, 仍 十 分 微小 . 表 1 中 数据 
的 负 号 , 代表 与 预 应 变 的 方向 相反 . 当 & 达 到 5% 时 ， 


| 


180| 后 
:\ 
sof SA 后 =58.41exp(-s/0.02832)+127.93 
We 
sse 2 
140[ ND 
ES 
万 =106.2008875 ~ a 
ee 
120 上 
0.00 0.05 0.10 0.15 


2 滞后 回 弹 瞬 时 模 量 及- 预 应 变 s 曲 线 


Fig.2 E, varied with pre-strain gi under strain e=1.7x10™ 


Ee 


弹 与 滞后 回 弹 


Tablel Initial springback and time-dependent springback without holding time 


Specimen oi/ MPa Ass 忆 /GPa Aé E:/ GPa Eu p/% 
EI-0.06 116 -0.00060 193.3 = a —0.00060 
PI-0.2 309 —0.00164 188.4 -0.00006 5150 -0.00170 3.53 
PI-5 593 —0.00434 136.6 -0.00020 2965 —0.00454 4.41 
PI-8 653 —0.00496 1317 —0.00028 2340 —0.00524 5.32 
PI-10 678 —0.00516 131.4 -0.00030 2260 —0.00546 5.49 
PI-15 729 -0.00578 126.1 —0.00042 1735 —0.00620 6.77 


Note: or 一 unloading stress, Aé:—initial springback strain, E.—initial modulus of time-dependent springback, 


Ags 一 time-dependent springback strain, E:— time-dependent springback modulus，su 一 total spring- 


back, p—time-dependent springback rate 


| 
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DP600 钢板 经 历 了 一 定 的 塑性 变形 后 , PI-5 试 样 的 
Ag 的 绝对 值 达到 0.02%, 比 PI-0.2 试 样 提 高 了 近 一 
个 数量 级 , 滞后 回 弹 现象 较为 明显 . 随 着 & 的 继续 提 
高 , A& 的 绝对 值 也 逐渐 增 大 . 当 & 为 15% 时 , DP600 
钢板 As 的 绝对 值 达到 最 高 值 0.042%, 比 PI-5 试 样 
提高 了 一 倍 多 . 定义 Ag 占 总 回 弹 (As+A 人 的 比例 为 
滞后 回 弹 占 比 P, 有 : 

DP=Ase/(As +As,) (5) 


图 3 为 DP600 钢 板 p 与 & 的 关系 . 由 图 3 可知 ,p 
随 着 & 的 提高 而 增 大 , p 由 PI-0.2 试 样 的 3.53% 增 加 
到 PI-15 试 样 的 6.77%, 其 变化 规律 接近 线性 分 布 . 
对 其 进行 线性 拟 合 为 : 

p=(21.927e, + 3.4287) x 10™ (6) 

图 4 为 PI-5, PI-8, PI-10 和 PI-15 试 样 钻 载 后 15h 
内 的 As 与 时 间 的 关系 曲线 . 由 图 4 可 知 , 试 样 的 A& 
随时 间 增 长 呈现 指数 衰减 . 在 滞后 回 弹 初 期 , 应 变 
可 弹 速 率 较 高 , 随 着 时 间 的 推移 , 应 变 回 弹 速 率 逐 
渐 降 低 . PI-5 试 样 滞后 回 弹 的 弛 和 豫 饱 和 时 间 约 为 
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4 滞后 回 弹 应 变 As 与 时 间 的 关系 曲线 


Fig.4 Measured time-dependent springback strain As var- 


ied with time 
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15000 s. 弛 豫 饱 和 时 间 随 着 有 的 增加 而 逐渐 增加 , 各 
曲线 基本 达到 回 弹 饱 和 时 , 饱和 弛 了 驳 时 间 约 为 
30000 s. 
4 室温 蠕 变 与 蠕 变 回 弹 

金属 室温 下 的 蠕 变 响应 取决 于 金属 内 部 组 织 
成 分 与 外 部 载荷 条 件 . 室温 蠕 变 可 以 理解 为 时 间 依 
赖 性 的 塑性 行为 , 并 且 通 常 在 低 于 宏观 屈服 应 力 
Ro 下 即 可 进行 “2 DP600 钢 板 室温 蠕 变 实 验 的 目 
的 有 2 个 : 一 是 获得 不 同 应 力 水 平 下 材料 的 室温 是 
变 规律 , 以 获得 黏 弹性 模型 中 的 柔 量 函数 构建 预 
测 模 型 趾 ; 二 是 获得 室温 里 变 回 弹 随 预 应 变 的 变化 
规律 . 
为 获得 蠕 变 柔 量 函 数 , 在 较 低 应 力 下 进行 了 不 
同 恒 应 力 蠕 变 实验 . 恒 应 力 分 别 为 309, 362, 371 和 
384 MPa, 其 在 图 1 中 的 真 应 力 -应 变 曲 线 上 , 分 别 对 
应 应 变 0.20%, 0.29%, 0.35% 和 0.41%. DP600 钢板 
在 上 述 较 低 应 力 下 保 压 4h 的 室温 蠕 变 曲线 如 图 $a 
所 示 . 可 见 , 随 着 恒 应 力 水 平 的 提高 , 蠕 变 应 变 逐 渐 
升 高 . 各 条 曲线 在 流 变 初 期 即 蠕 变 第 一 阶段 , 应 变 增 
长 较 快 , 应 变 率 较 高 , 之 后 随时 间 的 延长 逐渐 降低 , 在 
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保 压 4h 结 束 时 应 变 率 依次 为 0,2.7x10”, 5.1x10” 和 
1.4x108 s .可见 , 309 MPa 下 的 蠕 变 应 变 在 4h 内 基 
本 达到 饱和 , 这 与 蠕 变 柔 量 函 数 的 变化 规律 比较 符 
合 . 同时 , 蠕 变 应 变 的 数量 级 与 DP600 灌 后 回 弹 应 
变 均 为 10 因而 选取 309 MPa 下 的 蠕 变 曲 线 进 行 


义 Age 占 总 回 弹 (Asse+Aso) 的 比例 为 蠕 变 回 弹 占 比 v 
PC-0.3, PC-0.33, PC-0.4, PC-5, PC-10 和 PC-15 试 样 
的 .与 滞后 回 弹 实 验 的 类似, 随 着 应 力 的 提高 ， 
逐渐 降低 ; 而 EE 与 灌 后 回 弹 实验 的 EE 相 比 , 随 着 钊 
载 应 力 的 提高 呈现 良好 的 稳定 性 ; 而 vy 与 p 相 比 , 随 


拟 合 , 以 获得 蠕 变 柔 量 . DP600 钢 板 在 较 高 应 力 下 保 
压 4h 的 室温 蠕 变 曲线 如 图 5b 所 示 . 恒 应 力 选取 图 1 
中 应 变 5%, 10% 和 15% 所 对 应 的 真 应 力 593, 678 和 
729 MPa. 可 见 , 随 着 恒 应 力 水 平 的 提高 , 蠕 变 应 变 
也 逐渐 升 高 , 并 且 明 显 高 于 低 应 力 水 平 下 的 蠕 变 应 
变 . 在 较 高 应 力 水 平 下 蠕 变 实验 保 压 4h 结 束 时 的 
应 变 率 依次 为 9.7x10, 5.8x10” 和 1.6x10” s7. 

图 6a 为 DP600 钢 板 在 不 同 恒 应 力 下 的 蠕 变 划 
弹 曲 线 , 并 将 相关 实验 结果 列 于 表 2 中 . 可 见 , 随 着 
变 恒 应 力 Gi 的 逐渐 增 大 , 相同 回 弹 时 间 内 的 里 变 
时 回 弹 应 变 As. 和 里 变 回 弹 应 变 Ag. 的 绝对 值 也 
逐渐 增 大 . 定义 Gi. 与 Ass. 的 比值 为 蜂 变 回 弹 瞬 时 模 
量 E,, 其 物理 含义 为 单位 瞬时 回 弹 应 变 所 对 应 的 和 
载 应 力 ; 而 qi 与 As. 的 比值 为 蠕 变 回 弹 模 量 E, 其 物 


下 


着 伙 载 应 力 的 提高 略 有 降低 , 整体 变化 较为 稳定 ， 
因而 EE 和 v 均 呈现 了 一 定 的 稳定 性 . 由 图 6a 中 的 
PI-5, PI-10 和 PI-15 以 及 图 4 中 的 PC-5, PC-10 和 PC- 
15 试 样 的 实验 曲线 比较 可 知 , 经 过 4 h 的 保 压 过 程 ， 
Ag 的 绝对 值 与 As 的 绝对 值 相 比 明显 降低 . 将 PC- 
5, PC-10 和 PC-15 试 样 的 As 以 PL-$, PI-10 和 PI-15 
试 样 的 As 为 标准 , 定义 蠕 变 回 弹 降 幅 占 比 大 为 Ass 
的 减少 量 与 As 的 比值 : 

k=(Ae,— Ae,)/Ae, (7) 


k 与 & 的 关系 如 图 6b 所 示 . 可 见 ,k 与 & 旦 现 良 好 
的 线性 关系 , 线性 拟 合 结果 为 : 
上 = 0.143 +3.19e， (8) 


2.5 蠕 变 柔 量 


理 含义 为 单位 蠕 变 回 弹 应 变 所 对 应 的 外 载 应 力 ; 定 根据 顾 弹 性 理论 , 蠕 变 柔 量 J(D) 为 材料 在 单位 
(a) mn PC-0.2 ®。 PC-5 QF (b) 
+ 0.0F 各 v PC-0.4 。® PC-10 四 
9 ee A PC-15 可 
= pent We fe0.143+3.196 
S05 2 " 
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6 DP600 钢 板 蠕 变 保 压 4h 的 里 变 回 弹 应 变 和 里 变 回 弹 降幅 占 比 k 
Fig.6 Creep springback strain of DP600 steel sheet (a) and loss percentage k for creep springback (b) 
表 2 蠕 变 保 压 4h 后 的 蠕 变 瞬时 回 弹 与 蠕 变 回 弹 
Table 2 Initial springback for creep and creep springback after holding for 4h 
Specimen Goi/ MPa A ése EE./ GPa Aa。 E./ GPa Eu v/% 
PC-0.2 309 —0.00178 167.4 —0.00006 S5150 —0.00184 3.26 
PC-0.3 362 —0.00216 166.6 —0.00008 4525 —0.00224 3.57 
PC-0.33 371 —0.00224 165.1 —0.00008 4637 —0.00232 3.44 
PC-0.4 384 —0.00232 163.7 —0.00008 4800 —0.00240 3.33 
PC-5 593 —0.00446 132.9 —0.00014 4235 —0.00460 3.04 
PC-10 678 —0.00544 124.6 —0.00016 4237 —0.00529 必 :83 
PC-15 729 —0.00618 117.9 —0.00018 4050 —0.00636 2.83 


Note: ae 一 holding stress, Aé..—initial Springback strain for creep, E.—initial modulus of creep springback, 


As 一 creep springback strain, EF.—creep springback modulus, v—creep springback rate 
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阶 跃 应 力 下 的 应 变 响应 : 
J()=e(D)/oo (t=0) (9) 
式 中 , o% 为 所 加 恒 应 力 , s(0 为 材料 的 应 变 响应 . 由 于 
在 实验 中 不 易 实 现 瞬 态 加 载 , 不 同 加 载 条 件 下 的 里 
变 实验 均 具有 加 载 历史 , 加 载 历史 将 影响 黏 性 材料 
后 续 的 应 变 变化 25. 为 获得 不 含有 加 载 历史 效 应 
的 蠕 变 柔 量 , 同时 考虑 到 PC-0.2 试 样 的 恒 应 力 率 加 
载 , 采用 Prony 级 数 的 柔 量 拟 合 方法 进行 处 理 关 2 
一 维 黏 弹性 积分 型 本 构 方程 如 下 5: 
esO= oO+ [I -D600 (10) 
式 中 , 5 为 积分 变量 ; J(1- 人 为 经 过 时 移 和 反 转 变换 的 
蜂 变 柔 量 函数 ; c 为 应 力 率 函数 , 即 应 力 对 时 间 的 一 
阶 导数 . 注意 到 实验 初始 条 件 rc=0, 代入 到 上 式 有 : 
se()= [7JC-0od (11) 
以 外 载 时 刻 t 作 为 式 (11) 积 分 限 的 分 界 点 , 式 
(11) 可 改写 为 : 


sO=h -Dod + PI-Dod 0<n<D (12) 
注意 到 : 
1 7-9adi=0 (13) 
式 (12) 可 简化 为 : 
es0= 人 IJC-O6d (14) 


对 于 PC-0.2 为 恒 应 力 率 加 载 试 样 ，c 为 常数 ， 
式 (14) 可 进一步 写 为 : 
<sO=(o[ -0) de (15) 
为 了 便于 分 析 清 后 回 弹 , 选取 无 冲击 响应 的 4 
参数 Kelvin 链 进行 分 析 忆 , 柔 量 函 数 可 表示 成 : 
J()=(1 ~ exp(-tr N/E +( ~- exp(-t/r,)/E, (16) 
式 中 , 7 和 7, 分别 代表 第 一 个 和 第 二 个 Kelvin 单元 
的 延 滞 时 间 , 已 和 求 分 别 代 表 第 一 个 和 第 二 个 Kel- 
vin 单 元 中 弹簧 的 弹性 模 量 . 因此 , 在 阶 跃 应 力 下 的 
应 变 时 间 函 数 (不 含 瞬 态 响应 ) 为 : 
e()=00)()= 00( -exp(—t/T)NE + 
ol ~ exp(—1/7,))/E, (17) 
将 式 (16) 代 入 式 (15), 积分 得 计 入 恒 应 力 率 加 载 
历史 影响 的 应 变 函 数 : 
el(t) =(o /EL +(7T1/11)(1 一 exp(~t1/71))exp(-1/7)| 十 
(Co/E)TLE+(C2/)GL-exp(-n/ra))exp(-ur (18) 
由 于 PC-0.2 试 样 的 4=20 s, 将 其 代入 式 (18). 为 
简化 拟 合 , 采用 下 式 : 
a()=A, exp((-t+20)/7,)—A, exp(—t/7))+ 
A, exp((-t+20)/7,)-A, exp(-t/7,)+A, (19) 


式 中 ,4 和 4: 为 指数 项 的 待定 系数 , 4 为 常数 项 的 
待定 系数 . 采用 式 (19) 对 PC-0.2 试 样 的 蠕 变 曲线 
进行 拟 合 , 拟 合 结果 为 4o=0.00029, 4,=-0.00105， 
A:=-0.00826, 7,=107 以 及 7, =2283, 如 图 7 所 示 . 


对 比 式 (18) 和 (19), 可 得 : 
Ao=0o0/E,+oo/k, (20) 
Oo/E =A, Xi /7, QL) 
OolE,=A, Xt/7, (22) 


将 相关 系数 代入 式 (17), 可 得 DP600 钢板 在 阶 
跃 应 力 309 MPa 下 的 应 变 响 应 为 : 
el(l)= -1.9536x 10“ exp(-1/107) -0.7242x 10* x 
exp(-1/2283)+ 2.94559 x 10™ (23) 
对 应 的 蠕 变 柔 量 函数 为 : 
J()= -6.322x107exp(-U107)-2.3436x107x 
exp(—1/2283)+9.5338x10” (24) 
2.6 应 力 率 函数 
采用 中 心 差 分 数值 计算 方法 对 实验 数据 进行 
处 理 , 得 到 的 应 力 率 -时 间 关 系 曲 线 如 图 8a 所 示 , 对 
应 的 应 力 - 时 间 曲 线 如 图 8b 所 示 . 由 图 8a 可 知 , 应 
力 率 与 时 间 不 再 满足 线性 关系 , 所 以 加 载 历 史 效应 
&(D) 对 消 后 回 弹 的 影响 应 由 应 力 率 函数 与 柔 量 函 数 
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到 7 DP600 钢 板 柔 量 函 数 拟 合 及 局 部 放大 图 
Fig.7 Fitting curve of creep compliance for DP600 (a) and 


detail view of marked area in Fig.7a (b) 
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图 8 应 力 率 -时 间 曲 线 和 应 力 -时 间 曲线 


Fig.8 Stress rate-time curve (a) and stress-time curve (b) 


(ta is the starting point of stress rate function; ts, tc, 
to and ts are the times of unloading for specimens PI- 


5, PI-8, PI-10 and PI-15) 


卷 积 而 得 , 即 式 (14). 由 EI-0.06 试 样 的 实验 结果 可 
知 , 在 弹性 加 载 阶 段 , 室温 下 DP600 钢板 的 黏 性 行 
为 十 分 微弱 , 基本 无 滞后 回 弹 过 程 . 所 以 依据 图 8b 
所 示 , 将 计 入 历史 效应 的 应 力 率 起 始 时 刻 取 为 加 载 
后 第 40 s, 如 图 8a 和 b 中 点 所 示 , 对 应 应 变 约 为 
0.2%. 对 式 (14) 的 积分 下 限 修 正 为 a 


a0= [NO 40) 05) 


PI-5, PI-8, PI-10 和 PI-15 试 样 滞后 回 弹 实验 炙 
载 时 刻 分 别 为 , 如 各 和 也 如 图 8a 所 示 . 因而 , 对 于 
PI-5, PI-8, PI-10 和 PI-15 试 样 滞后 回 弹 实验 , 影响 滞 
后 回 弹 过 程 的 应 力 率 时 间 曲 线段 分 别 为 - 欠 段 t-te 
段 .t-tb 段 和 tw-ts 段 . 
2.7 不 同 预 应 变 滞后 回 弹 预测 

依据 式 (25), 通过 Matlab 数值 分 析 软 件 计算 应 
力 率 函 数 与 柔 量 函 数 离散 序列 的 卷 积 和 . 通过 将 卷 
积 和 拟 合 , 获得 PI-5, PI-8, PI-10 和 PI-15 试 样 的 加 载 
历史 效应 曲线 a(D), 如 图 9a 所 示 , 前 1000 s 的 放大 图 
如 图 9b 所 示 , 并 在 图 9(b) 中 标 出 了 PI-5, PI-8, PI-10 
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图 9 不 同 预 变形 量 的 卷 积 拟 合 曲线 及 局 部 放大 


Fig.9 Convolutions of historical loading effect under dif- 


ferent pre-strains (a) and detail view of marked area 


in Fig.9a (b) 


和 PI-15 试 样 对 应 的 扼 载 时 刻 , tc, 如 和 左 . 

由 以 上 讨论 可 知 ,在 PI-5, PI-8, PI-10 和 PI-15 
试 样 的 实验 中 分 别 为 ,tc,b 生 . 根据 式 (9), 秋 载 脉 
冲 效应 函数 -~&(D 可 表示 为 : 

-euO=-OiyO (26) 
式 中 , Ga 为 滞后 回 弹 实验 的 卸 载 应 力 , 式 (26) 中 的 负 
号 笈 不 印 载 过 程 的 应 疹 量 为 负 . 由 Boltzmann 著 
加 原理 , 分 别 在 图 9%b 中 所 示 的 加 ,tc 如 和 左 时 刻 , 时 
加 仓 载 脉冲 效应 函数 -es(D, 那么 弹性 后 效应 变 函 数 
&( 四 可 表示 为 : 


| 


Es(D)= et+t) -2,() (27) 
式 中 , &(1+) 为 经 过 时 移 变换 的 加 载 历 史 效 应 函数 . 
将 式 (25) 和 (26) 代 入 式 (27)， 


se,0=] + -人 一生 oN) 8) 


将 图 9a 中 的 卷 积 函数 和 柔 量 函 数 式 (24) 分 别 代 
入 式 (28), 可 得 PI-5, PI-8, PI-10 和 PI-15 试 样 的 弹性 
后 效应 变 曲 线 如 图 10a 所 示 . 为 便于 比较 不 同 试 样 
的 滞后 回 弹 应 变 , 将 弹性 后 效应 变 a(D) 减 去 零 时 刻 
的 初始 应 变 as(0) 可 得 滞后 回 弹 应 变 As(0O 的 预测 公 
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式 (29) 所 获得 的 各 曲线 清 后 回 弹 应 变 波形 如 
10b 所 示 . 可 知 , 式 (29) 较 好 地 反映 了 图 4 中 浪 后 
回 弹 应 变 的 绝对 值 随 着 a 增加 而 递增 的 趋势 . 图 11 
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图 10 弹性 后 效应 变 曲线 和 滞后 回 弹 应 变 预测 曲线 


Fig.10 Elastic after effect strain curves (a) and time-depen- 


dent Springback strain curves for prediction (b) 
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11 滞后 回 弹 实测 应 变 与 预测 应 变 随 预 应 变 s 的 变化 
规律 


Fig.1 Experimental strain and predicted strain of time-de- 


pendent springback varied with pre-strain & 
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式 , 有 : 为 滞后 回 弹 应 变 预 测 值 与 实验 值 的 比较 . 可 知 , 清 


后 回 弹 应 变 预测 值 较为 接近 实验 测量 值 , 并 且 随 着 
8 的 增加 均 呈 现 近似 线性 的 变化 规律 . 各 预测 曲线 
基本 达到 弛 豫 饱 和 时 , 饱和 时 间 约 为 500 s, 如 图 10b 
所 示 , 约 为 图 4 所 示 曲 线 的 1/60, 因而 弛 和 豫 饱 和 时 间 
的 预测 可 由 比值 160 近似 估 测 . 
考虑 到 蠕 变 保 压 过 程 对 蠕 变 回 弹 的 影响 , 由 式 
(8) 可 引入 因子 (1- 昌 , 与 式 (29) 相 乘 以 使 模型 计 入 保 
压 过 程 的 影响 , 有 : 
As D=As (DLL -月 (30) 
式 中 , As'Q) 为 计 入 保 压 过 程 的 蠕 变 回 弹 应 变 . 将 式 
(8) 代 入 式 (30)， 
As DO=As (D)(0.857-3.19s |) (31) 


| 


结论 
(1) DP600 钢板 在 滞后 回 弹 实验 中 , 滞后 回 弹 
十 模 量 随 着 预 应 变 的 增加 而 呈现 递减 规律 . 滞后 
弹 初 期 具有 较 高 的 回 弹 驱动 力 , 应 变 率 较 高 ; 随 着 
寸 间 的 推移 , 回 弹 驱 动力 减 小 , 应 变 率 逐渐 降低 ; 随 
预 应 变 的 增加 , 滞后 回 弹 占 总 回 弹 比 例 逐 渐 增 加 . 
(2) DP600 钢板 在 蠕 变 实验 中 , 蠕 变 回 弹 瞬 时 模 
随 着 预 应 变 的 增加 也 呈现 递减 规律 . 在 相同 的 
时 间 内 , 蠕 变 回 弹 应 变 比 直接 卸载 的 清 后 回 弹 应 
明显 减 小 ; 经 过 4 的 蠕 变 保 压 过 程 , 蠕 变 回 弹 模 
随 预 应 变 的 增 大 变化 较 小 , 接近 线性 关系 , 而 湾 
回 弹 模 量 的 变化 规律 呈现 减 小 趋势 . 

(3) 通过 对 线 黏 弹性 本 构 方程 积分 下 限 的 修正 ， 
构建 了 反映 非 线性 效应 的 新 型 滞后 回 弹 预 测 模型 . 
滞后 回 弹 的 预测 方向 与 实测 方向 相同 , 同时 滞后 匠 
弹 应 变 的 预测 值 接近 实测 值 ; 滞后 回 弹 应 变 预测 值 
的 绝对 值 随 预 应 变 的 增加 而 单调 递增 , 呈 线 性 分 
布 , 较 好 地 反映 了 滞后 回 弹 实 测 应 变 的 变化 规律 . 
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